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Abstract 

Unified theories containing a U(1)b-l gauge symmetry predict heavy Ma- 
jorana right-handed neutrinos. In such theories, cosmic strings may form at 
the B — L breaking scale. If the Higgs field forming the strings is also the 
Higgs field which gives mass to the right-handed neutrinos, there are right- 
handed neutrinos trapped as transverse zero modes in the core of the strings. 
When cosmic string loops decay, they release these neutrinos. This is an out- 
of-equilibrium process. The released neutrinos acquire heavy Majorana mass 
and decay into massless leptons and electroweak Higgs particles to produce a 
lepton asymmetry, which is converted into a baryon asymmetry via sphaleron 
transitions. 

PACS Numbers : 98.80.Cq, 12.10.-g, 11.30.Fs 
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Grand unified theories (GUT's) provide the standard scenario for baryogenesis, via the 
out-of-equilibrium decays of heavy gauge and Higgs bosons which violate baryon number 
(B) and lepton number (L). But it has been realized a decade ago that, unless the universe 
started with a non- vanishing B — L asymmetry, any B or L asymmetry generated at the 
GUT scale would be erased by sphaleron transitions An initial B — L asymmetry can be 
obtained in theories containing an extra gauge U(l) b-l symmetry via the out-of-equilibrium 
decays of heavy Majorana right-handed neutrinos This mechanism however requires 

either very heavy neutrinos or extreme fine tuning of the parameters in the neutrinos mass 
matrix ||. Also, the masses of the new gauge bosons must be bigger than the smallest heavy 
neutrino mass || . Hence there is a wide range of parameters for which the mechanism does 
not produce enough baryon asymmetry. 

In this letter, we show that in unified models involving an extra gauge U(l) b-l sym- 
metry, a primordial B — L asymmetry can be generated by the out-of-equilibrium decays 
of heavy Majorana right-handed neutrinos released by collapsing cosmic string loops. As 
a consequence of U(l) b-l breaking, cosmic strings may form at the B — L breaking scale, 
according to the Kibble mechanism ||. We call them B — L cosmic strings. The Higgs field 
mediating the breaking of B — L is the Higgs field forming the strings and it is the same 
Higgs field that gives heavy Majorana mass to the right-handed neutrinos. Hence, due to the 
winding of the Higgs field around the string, we expect right-handed neutrino zero modes 
@ trapped in the core of the strings. These zero modes are predicted by an index theorem 
||. There are also modes of higher energy bounded to the strings. We shall consider only 
the zero modes, which are the most favourable to be trapped. B — L cosmic string loops 
lose their energy by emitting gravitational radiation and rapidly shrink to a point, releasing 
these right-handed neutrinos. This is an out-of-equilibrium process. Right-handed neutri- 
nos acquire heavy Majorana mass and decay into massless leptons and electroweak Higgs 
particles to produce a lepton asymmetry. This lepton asymmetry is converted into a baryon 
asymmetry via sphaleron transitions. 

Topological B — L cosmic strings form when a gauge group G D U(1)#-l breaks down 
to a subgroup H \J(1) b-l of G, if the vacuum manifold ^ is simply connected, that is if 
the first homotopy group is non-trivial. If = I but string solutions still exist, 

then embedded strings || form when G breaks down to H. Embedded strings are stable for 
a wide range of parameters. In left-right models, embedded B — L strings usually form. In 
the simple U(l) extension of the standard model SU(3) C x SU(2)i x U(l)y x U(l)y where 
Y' is a linear combination of Y and B — L topological strings form [|10j. In grand unified 
theories with rank greater than five, such as SO(10) or E(6), B — L cosmic strings may form, 
depending on the symmetry breaking pattern and on the set of Higgs fields used to do the 



breaking down to the standard model gauge group |TT],|T0[- There is a wide range of theories 
which contain both U(1)b_l and B — L cosmic strings. 

Consider a unified model with an extra U(l) b-l gauge symmetry and stable B — L cosmic 
strings. The gauge and Higgs fields forming the strings will be the B — L associated gauge 
boson A' and the Higgs field 4>b-l used to break \]{1) B -l- Right-handed neutrinos acquire 
heavy Majorana mass via Yukawa couplings to 4>b-l- The U(1)b-l part of the theory is 
described by the Lagrangian 

L = \f»„r + (D^B- L )KDy B -L) - V^b-l) (1) 
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+iV L ^D„v L + wrI^D^Vr (2) 

+ i\4>B-L^R^ C L -iM*B-L V L P R + L f ■ ( 3 ) 

The covariant derivative = — iectf' where e is the gauge coupling constant and Y' is 
a linear combination of B — L and Y . A is a Yukawa coupling constant, and V(4>b-l) is 
the Higgs potential. The spinor N = vr + v c L is a Majorana spinor satisfying the Majorana 
condition N c = Cj^N* = N, where C is the charge conjugation matrix. Hence N has only 
two independent components, two degrees of freedom. Lf is the fermionic B — L part of the 
Lagrangian which does not contain neutrino fields. 

For a straight infinite cosmic string lying along the z-axis, the Higgs field 4>b-l and the 
Y' gauge field a in polar coordinates (r, 9) have the form 

= f(r)e md (4) 

ag = —nr a z = a r = (J, (5) 

er 

where n is the winding number; it must be an integer. Most strings have winding number 
n — 1; strings with winding number \n\ > 1 are unstable, r is the string's generator; it 
is the normalised U(l)y generator. It has different eigenvalues for different fermion fields. 
The functions f(r) and g(r) must satisfy the following boundary conditions 

/(0) = and / -> t] B -l as r -> oo, (6) 
5^(0) = and g — * 1 as r — > oo , (7) 

where ??,b-l is the scale of B — L breaking. The exact forms of the functions f(r) and g(r) 
depend on the Higgs potential V^b-l)- 

From the Lagrangian (|3|) we derive the equation for the right-handed neutrino field: 

vfDyfi, - iX^B-itl = (8) 

where v c L = C^qV r . Solving (H), we find that Majorana neutrinos trapped as tranverse zero 
modes in the core of B — L cosmic strings have only one independent component. For an 
n — 1 vortex it takes the form : 



N 1 = P(r,9)a(z + t) (9) 

where /3(r, 6) is a function peaked at r = which exponentially vanishes outside the core of 
the string, so that the fermions effectively live on the strings. The z and t dependence of a 
shows that the neutrinos travel at the speed of light in the — z direction, so that they are 
effectively massless. In an n = -1 vortex, the function a = a(t — z), so that the fermions 
travel at the speed of light in the +z direction. These fermions can be described by an 
effective theory in 1 + 1 dimensions. The usual energy to momentum relation 

E = P (10) 

holds. We have no boundary conditions in the 1 spatial dimension, and the spectrum of 
states is continuous. In the ground state the Fermi momentum of the zero modes is pp = 0. 
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The field solution (|]) and the energy to momentum relation ([10]) have been derived for 
fermions on a straight infinite string. However physical cosmic strings are very wiggly and 
are not straight. Hence relations and ([TOD do not hold in the physical case. Neither do 
they hold for cosmic string loops, even if the latter are assumed to be smooth. On a cosmic 
string loop, fermions are characterised by their angular momentum L. The energy relation 
becomes [Ol 



E = ^±jl = P + ± (ii) 

R 2R K J 



and hence the energy spectrum is 



fn+ r 



E =—^ (12) 



where n G N. We see from Eqs. (|T0"|), (|TT|), and ([12]) that, when R is very large, the string 
looks locally like a straight string. We have an almost continuous spectrum of states. The 
fact that the string gets a finite curvature acts as a perturbation on the string bound states. 
The energy levels get quantised and the Fermi energy gets a non- vanishing value, Ep = 
As the string loop shrinks, its radius R decreases and we see from Eq.([T^) that the Fermi 
energy increases and that the separation between energy levels gets wider. 

Assuming that a loop decays when its radius R becomes comparable to its width w ~ 
Vb-li we deduce that the Fermi energy level when the loop decays is E F ~ \t\b-l-, where 
the B — L breaking scale t]b-l is of the order of the right-handed neutrino mass. Ep is lower 



than the energy needed by the neutrinos to escape the string [0. Hence, when a cosmic 
loop decays, it releases at least n u = 1 heavy Majorana neutrinos. Quantum fluctuations 
and finite temperature corrections may increase n v . Part of the final burst of energy released 
by the decaying cosmic string loop is converted into mass energy for the gauge and Higgs 
particles released by the string, and into mass energy for the neutrinos. A decaying B — L 
cosmic string loop releases heavy B — L Higgs particles which can decay into right-handed 
neutrino pairs, and hence increase n u . This is an out-of-equilibrium process. Due to angular 
momentum conservation, the massive Majorana neutrinos released by a decaying cosmic 
string loop which were trapped as transverse zero modes are spinning particles. 

Heavy Majorana right-handed neutrinos interact with the standard model leptons via 
the Yukawa couplings 



hijliH^i/Rj + h.c. (13) 



where I is the usual lepton doublet; for the first family / = (e, H ew is the standard model 
doublet of Higgs fields. Majorana right-handed neutrinos can decay via the diagrams shown 
in Fig. l.a. and l.b. CP is violated through the one loop radiative correction involving a 
Higgs particle as shown in figure l.b. The right-handed neutrinos are out-of-equilibrium, 
and hence a lepton asymmetry can be generated. The lepton asymmetry is characterised 
by the CP violation parameter e which, assuming that the neutrino Dirac masses fall into a 
hierarchical pattern qualitatively similar to that of the leptons and quarks, is estimated to 
be! 

TTV M N2 K 1 
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where 771d 3 is the Dirac mass of the third lepton generation, v is the vacuum expectation 
value of the electroweak Higgs field v =< H ew >= 174 GeV, M N1 and M N2 are the right- 
handed neutrino Majorana masses of the first and second generation respectively and 5 is 
the CP violating phase. 

The corresponding B—L asymmetry (we use B—L instead of L since the (S+L)-violating 
electroweak anomaly conserves B—L) must be calculated solving Boltzmann equations which 
take into account all B, L and B + L violating interactions and their inverse decay rates. We 
can however calculate the B — L asymmetry produced taking into account only the out-of- 
equilibrium decays of right-handed neutrinos released by decaying cosmic string loops and 
assuming that the rates of inverse decays are negligible. Hence an upper limit on the baryon 
number per commoving volume at temperature T is then given by (13 



B(T)^-Mh, (15) 

where s is the entropy at time t and N v (t) is the number density of right-handed neutrinos 
which have been released by decaying cosmic string loops at time t. Recall that the temper- 
ature T is related to the cosmic time t via the relation t = 0.3 g* 2 where g* counts the 
number of massless degrees of freedom in the corresponding phase and M p i is the Planck 
mass, s, the entropy at time t, is given by s = ^7r 2 ^*T 3 . N v (t) is approximately n v times 
the number density of cosmic string loops which have shrinked to a point at temperature 
T. Assuming that sphaleron transitions are not in thermal equilibrium below T ew , and ne- 
glecting any baryon number violating processes which might have occured below T ew , the 
baryon number of the universe at temperature T < T ew is then given by 

B = B(T ew ), (16) 

which is also the baryon number of the universe today. If sphaleron transitions are also 
rapid below T ew , we should include the neutrinos released below T ew . However, below T ew 
the number density of decaying cosmic string loops is negligible, and hence this would not 
affect the result in any sense. 

The number density of decaying cosmic string loops can be estimated from the three 
scales model of reference [15]]. The model is based on the assumption that the cosmic string 



network evolution is characterised by three length scales £(£), and x(t) related to the 
long string density, the persistence length along the long strings (which is related to the fact 
that the typical loop size is much smaller than £), and the small scale structure along the 
strings respectively. 

Cosmic string loops lose their energy by emitting gravitational radiation at a rate E = 
—FioopsGfi 2 JL5[], where fi ~ T c 2 is the string mass-per-unit-length and T c = t]b-l is the 
critical temperature of the phase transition leading to the string network formation. G is 
Newton's constant. The numerical factor T [oops ~ 50 — 100 depends on the loop's shape and 
trajectory, but is independent of its length. The mean size of a loop born at £& is assumed 
to be (k — l)Ti oops Giitb- At a later time t, it is then Ti oops Gfi(ktb — t). The loop finally 
disappears at a time t = kty,. Numerically, k is found to lie between 2 and 10. 

The rate at which the string loops form in a volume V is given by [|T^ 

^ ^ (k — l)Ti oops Gfitf ^ ^ 
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where the parameter v can be expressed in terms of the various length scales of the model 
which vary with time. We start with £ ~ £ ~ \. Then £ and later £ will start to grow 
and will evolve to the scaling regime characterised by £(£) and £(£) ~ t. The length scale \ 
grows much less rapidly. Therefore v varies with time. In the scaling regime, v is estimated 



to lie in the range v = 0.1 — 10 ||15|| . 



Note that it has recently been shown that cosmic string networks reach the scaling 
solution at a time t* much smaller than previously estimated |l4j|. The authors of ref. WA 



find that in the radiation dominated, for minimal GUT strings, t* = 8 x 10 2 £ c [H|, where 
t c is the time at which the strings form, and the associated temperature is T* ~ 10 14,5 GeV. 
Hence our approximation for the rate of cosmic string loop formation is suitable; it leads to 
a lower bound on the number of cosmic string loops. Only numerical simulations could lead 
us to a better estimate, but it is beyond the scope of this paper. 

Since it has been shown that most of the loops formed have relatively small size, we 
shall assume that the number of loops rejoining the network is negligible and thus that the 
number of decaying loops is equal to the number of forming loops. Hence the number density 
of right-handed neutrinos which have been released by decaying cosmic string loops at time 
t is given by 

where where r(t) is the cosmic scale factor and n v is the mean number of right-handed 
neutrinos released by a single decaying loop. In the radiation dominated era r{t) ~ t?. 
After integration we obtain 

6 [k - L)L loopsG/j, (0.3)3^2 M pl M pi J 

Hence the baryon number per commoving volume today produced by the decays of heavy 
Majorana right-handed neutrinos released by decaying cosmic string loops given in Eq.fllED 
becomes 

(0.3 7^ (k-l)klT loops M pl M N2 (2Uj 

where we have used Eq.(p^) for the CP violation parameter e. The produced B asymmetry 
depends on the cosmic string scenario parameters, on the neutrino mass matrix and on the 
strength of CP violation. 

We now calculate the lower and upper bounds on B which correspond to different values 
of the parameters in the cosmic string scenario. We fix the neutrino mass matrix parameters 
and assume maximum CP violation, i.e. sin 5 = 1. We assume that the Dirac mass of the 
third generation fermions lies in the range m£> 3 = 1 — 100 GeV, and that the ratio of the 
right-handed neutrino masses of the first and second generation -jg^ = 0.1. With the above 
assumptions, we find B to lie in the range 

B ~ (1 x 10- 10 - 5 x 10- 2 ) gl (-^), (21) 
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and we see that B strongly depends on the cosmic string scenario parameters. Recall that 
the baryon number-to-photon ratio ^ is related to the baryon number of the universe B by 
jf- = 1.80 g*B. Hence, if g* ~ 10 2 , our mechanism alone can explain the baryon-number- 
to-photon ratio predicted by nucleosynthesis, ^ = (2 — 7) x 10~ 10 , with the B — L breaking 
scale in the range 

Vb _ l = (2 x 10 7 - 3 x 10 16 ) GeV. (22) 

We point out that the result could be better calculated solving Boltzmann equations, which 
take into account all B, L, and B + L violating interactions and do not neglect the inverse 
decay rates. Furthermore, the rate of decaying cosmic string loops can be calculated via 
numerical simulations, which would have to take into account the different regimes of the 
network evolution which occur during and after the friction dominated era. This may change 
the allowed value for the B — L breaking scale by a few orders of magnitude. Note also that 
if CP is not maximally violated, the B — L breaking scale will be shifted towards higher 
values. Finally, we recall that when a cosmic string loop decays, it also releases massive Higgs 
bosons (pB-L and massive gauge bosons a which can decay into right-handed neutrinos. This 
process is not taken into account here because the masses of the Higgs and gauge bosons 
and of the neutrinos are very close to each other and the Higgs and gauge fields can also 
decay into other particles. 
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FIG. la. The tree-level diagram for 
right-handed neutrino 
decay. 



FIG. lb. The one loop diagram for 

right-handed neutrino 
decay. 
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